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We present a comprehensive full-sky 3-dimensional analysis of the weak-lensing fields and their 
corresponding power spectra. Using the formalism of spin-weight spherical harmonics and spherical 
Bessel functions, we relate the two-point statistics of the harmonic expansion coefficients of the weak 
lensing shear and convergence to the power spectrum of the matter density perturbations, and derive 
small-angle limits. Such a study is relevant in view of the next generation of large-scale weak lensing 
surveys which will provide distance information about the sources through photometric redshifts. 
This opens up the possibility of accurate cosmological parameter estimation via weak lensing, with 
an emphasis on the equation of state of dark energy. 

PACS numbers: 


I. INTRODUCTION 

Like the Cosmic Microwave Background (CMB) a decade ago, the study of gravitational weak lensing on a cosmic 
scale is now entering a promised golden age. As with the CMB, it is well understood theoretically it has now 

been detected (see e.g. and ambitious experimental projects are under way. The uniqueness and thus 

the appeal of weak lensing lies in its clean and direct sensitivity to the total mass distribution of the Universe. Indeed, 
it arises from the deflection of light by the gravitational potential of the matter along the photon path, regardless of 
the precise nature or state of the intervening matter. With weak lensing one can hope to measure for the first time the 
3D matter power spectrum over a wide range of scales, in principle independently of any model for evolution. It thus 
provides complementary information to that rendered by the CMB alone or other lower redshift cosmological probes. 
The CMB does provide the cleanest possible window into the physics of the early Universe and its parameters iO, 
but lacks strong sensitivity to the Universe’s subsequent evolution. Large-scale galaxy surveys are hindered by the 
presence of non-linearities which affect the luminous baryonic matter on small scales and have to rely on assumptions 
about (or deductions of) the relationship between luminous matter and mass distribution on large scales flit (ill fl3l] . 
Supernova studies are afflicted by complicated systematics Hi- Hence, weak lensing can considerably help in the 
determination of poorly constrained parameters linked to the matter evolution and the moderate-redshift evolution of 
the redsliift-distance relation. The most significant parameter is the equation of state of dark energy, w = p/pc 1 , and 
its time evolution. As is now generally accepted, there is strong evidence from a variety of cosmological probes like 
supernovae [l5| and the CMB [kJ that a dark energy component makes up about 70 % of the total energy density 
budget of our Universe but its precise nature and state remain a mystery. Answers are unlikely to come from particle 
physics and our hopes lie in cosmology. 

A variety of deep, large sky weak lensing surveys are either ongoing, scheduled or planned in the foreseeable future 
(Deep Lens Survey fH}, NOAO Deep survey |l7|, The CFHT Legacy survey |Ts|, darkCAM on VISTA, The Panoramic 
Survey Telescope and Rapid response System (Pan-STARRS) 20j, the Large Synoptic Survey Telescope ;21j and 
SNAP 0|)- ^ one can rigorously account for the various systematics (such as the seeing, calibration, anisotropic 
point spread functions, redshift uncertainties etc.) likely to permeate the data, such surveys can offer if not definite 
answers at least fundamental clues on the nature of the dark energy, which lensing surveys are beginning to constrain 
w |2^| . In combination with other data sets, weak lensing surveys data can further overcome intrinsic degeneracies 
characteristic of each of the cosmological probes and strengthen the constraints on our cosmological model. In 
particular, predictive studies have been done which explore the possible improvements for parameter estimation of 
considering the three dimensionality of weak lensing, where photometric redshifts are used to provide estimates in three 
dimensions of the weak lensing shear field. As shown by Ishak 12 .'ll . cosmic complementarity and 3D lensing tomography 
bring uncertainties on the equation of state to the level of a few percent. The use of full distance information on 
weak lensing surveys, enabling a full 3D analysis rather than tomography, can further reduce statistical errors on 
cosmological parameters [24j . 

Besides reducing statistical errors, there are other reasons for wishing to have photometric redshifts: weak lensing 
studies have to deal with systematic errors which could dominate the error budget as surveys improve (see 0 and 
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references therein). The most important of these is the redshift distribution of the sources but correlated source 
galaxy intrinsic ellipticities (to date still not well-known) and the effect of source clustering are likely to become 
a considerable nuisance in the future. A reasonably accurate knowledge of the redshift of the sources can help in 
understanding, quantifying and ultimately removing such systematics. Furthermore, since acquisition of photometric 
redshifts requires several images of each galaxy in different bands, there is scope for better shape measurement, 
independent lensing studies, and so on. 

In the past, most of weak lensing analysis were essentially limited to 2D approaches due to the lack of, or imprecise, 
photometric redshift information about individual galaxies. Moreover, such studies were usually applied to small 
patches of the sky necessarily involving a flat-sky approximation (see e.g. |25|, ES and 130 for an application to 
data). A useful exception to this early trend is the theoretical work by Stebbins [23 who explored the 2D curved sky 
by means of the tensor spherical harmonics formalism. But over the last few years or so, with the arrival of surveys 
with more accurate photometric redshift information and the prospect of them reaching wider areas of the sky, 
interest has steadily grown in exploring the inherently full-sky 3D information contained in the lensing observables 
(see e.g. 0 0 So] ). The way of dealing with any source redshift information has thus evolved over the years 
from basically a black and white picture to a ‘fabulous Technicolour’ approach. Originally it was used as a way of 
determining the redshift distribution of the sources only. It then evolved to a so-called tomographic analysis, where 
the sources were divided up in slices at different redshifts and a 2D anal ysis was then performed on each one of the 
slices. This can be useful when trying to detect clusters of galaxies liul IdTl liiil [dill. As shown by Hu [34j at a statistical 
level such method also presents certain advantages depending on the parameters one is interested in retrieving. For 
instance, for the amplitude of the matter power spectrum there can be significant gains when a source population is 
split into two, but little is gained by further finer subdivisions. For w the gains are much larger as the 2D analysis 
constraints are weak. It is worth noting here that other 2D routes have been explored in pursuit of cosmologically 
sensitive weak-lensing estimators such as the study of higher-order statistics (for bispectrum detections s ee 3> ~). (36)3 
but if one wants to place strong constraints on re a 3D study is compulsory. Several theoretical studies |37l I38L I.'IPL 14(1) 
investigated how the lensing equations can be inverted to extract the 3D gravitational potential directly, thus offe ring; 
the possibility of reconstruction of the 3D mass density field from weak-lensing data. Only recently in Heavens |24l . 
was a truly 3D spectral statistical analysis on the full-sky developed, in which the individual source redshifts were 
taken into consideration per se. In there, the attention was mainly focused upon understanding the amount of extra 
information on the cosmological parameters that can be gleaned from a 3D shear map, rather than on establishing 
a framework for 3D weak lensing studies. Motivated by the previous theoretical considerations and the upcoming 
experimental prospects, it is timely to develop in detail a description of the weak lensing observables on the full three 
dimensional sky. 

A full-sky 3D study of the weak lensing observables needs to combine at its root both the 3D weak lensing char¬ 
acteristics and the full-sky formalism. The weak lensing effect can be observed as a local modification of the surface 
number density of background galaxies (the magnification) and a distortion of their shape (the shear). As has been 
pointed out repeatedly before, the weak shear components present striking similarities with the Stokes parameters Q 
and U describing the linear polarisation of the CMB light, the major difference being its additional radial dependence. 
Like them, they are not invariant under a rotation of the reference frame used to describe them but rather transform 
as rank-2 tensors. One can then legitimately extend to the third radial dimension past work done for the full-sky 
polarisation of the CMB 0, |43, |43j, 0] and adapt it to the weak lensing analysis. In this work, we develop in more 

detail Heavens’ original 3D full-sky analysis |241 weaving into our approach both the CMB polarisation formalism and 
past weak lensing theoretical 2D studies [23, 12a, 123 ■ We use a spectral decomposition as it allows one to restrict the 
analysis to, for example, the linear or mildly non-linear regime and rely on the spin-weighted formalism, originally 
developed by Newman and Penrose in the 1960s [4^ . 

This paper is organised as follows. In section un we review the main results of the theory of gravitational weak 
lensing which we generalise to 3 dimensions so that the relevant theoretical expressions depend explicitly on the radial 
distance from us. In subsection ITTbI we start by presenting the weak lensing theoretical expressions in the familiar 
tensorial formalism, both in the full-sky and in the flat-sky approximation. In subsection irm we introduce an 
alternative differential formalism we name the edth formalism, which encompasses the full and the flat-sky geometries. 
Such a mathematical apparatus reveals itself to be extremely convenient when working on the full-sky and has been 
extensively used in studies of the CMB polarisation but which we develop here for 3D weak-lensing fields. In section lTTO 
we spectrally decompose in the 3D full-sky the (spin-weight 2) weak lensing shear and (scalar) convergence fields in 
terms of spin-weighted spherical harmonics and Bessel functions. We derive their expansion coefficients in function 
of the gravitational potential coefficients by making full use of the edth formalism. In subsections IIV Al and IIV Bl we 
calculate the weak lensing 3D power spectra in the full-sky and derive small-angle limits for completeness. Finally in 
section El we present results for the 2-point correlation function of the 3D shear field. In the Appendixes, we review 
in more detail the notation, the mathematics and the conventions chosen of the spin-weight s functions and of the 
associated geometrical spin raising and lowering (3 and 3) operators defined over any two-dimensional Riemannian 
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manifold. 


II. WEAK LENSING THEORY IN 3D 


A. The lensing potential 


In the linear regime, many scalar fields on the sky that are associated with large-scale structures can be interpreted 


as line-of-sight integrations of functions of the gravitational potential $ with a given wi ' ’it. A few examples of such 
scalar fields are the Integrated Sachs-Wolfe effect [ZtJ or the Ostriker-Vishniac effect imprinting themselves on 
the CMB and the gravitational weak lensing by cosmological structure (for a comprehensive review see Q)- For the 



gravitational weak lensing case, one can associate the so-called weak lensing potential (f> for a given source at a 3D 


position in comoving space r = (r, 9, ip) to the peculiar gravitational potential $ defined along the line-of-sight via 


(see Q) 



(1) 


where the Born approximation was assumed (i.e. the path of the integration, corresponding to the path of the photons 
emitted by the source, is assumed to be unperturbed by the lens). Although usually the lensing potential cj) is regarded 


as a 2D radial projection of the 3D gravitational potential, it is in reality a 3D quantity. It is customary to average 
over the redshift distribution of the source galaxies, but this is not necessary if one has distance information about the 
individual sources. Here and in the remainder of the manuscript, bold letters denote 3-dinrensional vectors, c is the 
speed of light, r = r{t) is the comoving distance of the source at instant t from the observer at the origin (r = 0) and 
is the comoving transverse dimensionless separation for points separated by ch A We have /A(r) = sinr,r 
and sinhr for a closed ( k = 1), flat ( k = 0) and open ( k = —1) universes respectively. The gravitational potential *I> 
is related to the underlying overdensity field 5(r) = 5p(r)/~p by the Poisson equation 



( 2 ) 


where the 3D gradient V r is defined relative to comoving coordinates, is the present-day total matter density, Hq 


is the Hubble constant today in units of km/s/Mpc and a(t) = 1/(1 4- z) is the scale factor. We use the comoving 
gauge or total-matter gauge (see e.g. pfl l50 j). 


If one wants to make a spectral expansion, such as a Fourier transform, then immediately one has a subtlety to 
consider. The lensing potential field, indirectly observed through its effect on the shapes of lensed galaxies, is not 
homogeneous, as it is viewed on the past light cone of the galaxies, and is a function of the gravitational potential 
which evolves with cosmic epoch. When referring to the transform of a field, such as the gravitational potential, at a 
lensed galaxy at distance r, we will mean the transform of the homogeneous field existing everywhere at the cosmic 
epoch corresponding to the time the observed light left the distant galaxy. Thus the coefficients of the expansion 
depend on the lookback time of the observation, and hence, rather paradoxically, on the distance r itself. 

There is a natural choice of basis functions to use for an expansion. Spherical coordinates are natural for various 
reasons, partly because the lensing potential is a radial integral, partly because partial sky coverage is more easily dealt 
with, as demonstrated in large-scale structure studies by |5lj|. In addition, errors in distances (from using photometric 
redshifts as distance indicators, for instance) are radial errors. The first of these means the coefficients of the expansion 
of the gravitational potential and of the lensing potential are straightforwardly related. The choice of basis functions 
is motivated by Poisson’s equation ©. It makes sense to use the eigenfunctions of the Laplacian operator, since then 
the coefficients of the gravitational potential and of the density field are closely related (essentially by a factor —k~ 2 ). 
In Cartesian coordinates, the eigenfunctions are the familiar exponential functions of Fourier analysis. If however 
the Laplacian is written in spherical coordinates (r,9,ip), then in flat space the eigenfunctions become products of 
spherical harmonics and spherical Bessel functions: Ye m (0,ip)jc(kr), with eigenvalue — k 2 . So for scalar fields /(r) in 
a flat background geometry, the natural 3D expansion is 



( 3 ) 


where the numerical factor and the presence of k are chosen for convenience. The inverse transform is 



f( r ) = 


( 4 ) 
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This is readily obtained from the orthonormality of the spherical harmonics and the orthogonality of the spherical 
Bessel functions (see Eqs. HB2I) and m)- 

Using this 3-dimensional expansion, one can relate the transform of the lensing potential to the gravitational 
potential field in a flat geometry [24| by 






drrje(kr) 




je{k'r')^e m (k'\r'). 


(5) 


It is worth remarking that our expression differs from Eqs. ( 6 ) and (7) found in Heavens [2] due to a different choice 
of the 3D expansion conventions. One can further relate r) to the expansion of the matter overdensity r) 

by using the Poisson equation: 


$tm{k\ r) 


3 

2 k' 2 a(r) 


Se m (k;r), 


( 6 ) 


where here we put the time dependence (see discussion above) in the coefficients of the gravitational potential and 
matter overdensity fields by explicitly writing the coefficients as functions of the distance r in addition to £, m and 
k. We point out that the expansion coefficients of the lensing potential (f) do not have the time dependence explicitly 
shown because (f> is not by definition a homogeneous and isotropic field in 3D space. 

These equations establish the relations between the scalar fields <(>, d> and 6. Since the statistical properties of S are 
known for a given cosmological model, this opens up the possibility of using weak lensing in 3D (via <p) for estimation 
of cosmological parameters. This would normally be done using the weak lensing shear field, which is not a 

scalar field, but rather a spin-weight 2 held. We will consider the statistics of this observable held in Section III. 


B. Weak lensing and the gravitational potential in the tensorial formalism 

In this subsection we review the main results of the theory of gravitational weak lensing in the tensorial formalism 
on the 3 dimensional spherical sky. In addition, we introduce the relevant formulae in the 3 dimensional hat-sky 
approximation for completeness. 


1. Weak lensing on the full-sky 

The 2 dimensional distortion of images of distant sources, located at a certain 3D comoving position in space r, 
caused by the weak gravitational lensing by intervening structures is given by 0, |Z2, H2, E3 

[ v * V j - ^iV 2 ]<Kr) = [7i( r ) cr 3 + (7) 

where the indices i. j = (0, p ) stand for the polar 2D coordinate indices on the sphere, ay and er 3 are the Pauli matrices 
dehned on the 2D spherical sky, gij is the 2D metric on that surface given by g = diag( l,sin 2 0) and (j> is the lensing 
potential related to the gravitational potential <I> by Eq. © Here, and throughout, 71 and 72 are the components of 
the weak lensing shear produced by the gravitational tidal held, and which may be conveniently written as a complex 
shear y(r) = 71 (r) + ^(r). They correspond to the two orthogonal modes of the distortion which are, a priori, 
measurable on the sky with respect to a chosen hxed coordinate system. At hxed r, comparison of 7 at different 
locations ( 6 , p) on the sphere is not possible. In curved space, comparisons are only possible locally, and, as 7 is not a 
scalar, parallel-transport of 7 from place to place is path-dependent. This fact lies at the source of the interest in hrst 
developing rotationally invariant CMB polarisation components in terms of electric and magnetic components l4ll| 
as we will see shortly in the context of weak lensing. The differential operator V corresponds to a 2D covariant 
derivative on the sphere of radius r, which for a scalar is given by the partial derivative Vi = <9i, for a covariant 
vector Xi (i.e. rank-1 covariant tensor) is given by V jXi = djXi — T^X^ and for a contravariant vector X 1 is given 
by VjX 1 = djX 1 + T l kj X k with the Christoffel symbol PC depending on the metric g^ (see for instance [54|). In our 
case, for which g = diag( l,sin 2 0), we will have only three non-zero Christoffel symbols which are T 31 = — sin0cos0 
and Tj-j = T 30 = cot0. By construction, the Christoffel symbol is always symmetric in its lower indices T k - r = T-t. 
Due to this symmetry property, there is an important result for scalar fields X which is that V^V^-A = VfcV,;A. In 
the remaining of this work we will apply this property to the lensing potential <fi. 
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The lensing potential (f> can alternatively be observed by means of the magnification via the isotropic convergence 
scalar field n defined as 


[«( r )h = « {r)Iij 




( 8 ) 


where V 2 = V;V* and I{ 3 is the Identity matrix. We name [n]i 3 the convergence field tensor. Note that the weak 
lensing regime corresponds to lensing configurations such that |«| <C 1 and lyj <C 1. In this specific regime, the 
magnification and the distortion of galaxy sources are so small that one cannot measure them individually, rather one 
needs to perform a statistical study of the lensed population. 

It will be useful later to have the expressions 0 and 0 explicitly expanded in terms of covariant derivatives on 
the 2D spherical sky. Transforming from the Cartesian {e x (n), e y (n)} to a spherical polar 2-dimensional coordinate 
system {eg(n), e v (n)}, where we choose the basis of the two coordinate systems to be aligned such that dx = dd and 
dy = sin 9dtp, we can explicitly express Eq. 0 for the weak lensing shear as 


[7 (r) ]ij 


( 7 i ( r ) sin 6 » 7 2 (r) 

\sin# 72 (r) — sin 2 $71 (r) 


^[VeV,} - esc 2 V^Ve 

V v V e ^[V^V^ - sin 2 dWgWe] 


(j){r). 


(9) 


[ 7 \ij is the 2x2 symmetric and traceless cosmic shear tensor field associated to a source at a 3D position r defined 
on the sky and again Vi corresponds to a 2D covariant derivative on the sphere. We used the general result for scalar 
fields VeV v <f> = V v Vg The sin# terms in the weak lensing shear tensor [ 7 ]^ appear because the orthogonal (9, <p) 
basis is not orthonormal. This result has been repeatedly used when describing the polarisation tensor V, 3 of the 
CMB on the full-sky (see for example : 55j). One can extract from Eq. 0 the individual components of the weak 
lensing shear 71 and 72 


7 i( r ) = ^ [V$ Vg esc 2 #V^ V cj)(r ), 

72(f) = esc 9'VcpVg 4>{r). (10) 

Similarly, one obtains for the convergence field tensor [k] i 3 given by Eq. 0 the following expanded expression 


= ( l sin 2 0 ) K ^ = 


= 2 


+ esc 2 


0 

+ sin 2 0S7gS7e] 


4 >{ r ) 


(ii) 


where k, the scalar convergence field in the 3D full-sky, is given by 

k (r) = - [V e V s + esc 2 V v ] <j>(r). 


( 12 ) 


2. Weak lensing on the flat-sky 


When analysing data on a small patch of the sky, one can use the flat-sky, or small-angle, approximation. There, one 
defines the standard Cartesian coordinate system {ea,(n), e y (n)} with metric given by g = diag( 1,1). The covariant 
derivatives simplify to the standard partial differentiation operators Vi —> di. The above Eq. 0 then reduces to j 2 (| 

[7 (r) ]« = ( ) = (d t d 3 Is* V 2 ) <j)(r) (13) 

in such a way that we recover the well-known expressions for the components 71 and 72 of the shear on the 3D flat-sky 

7i( r ) = \{dl~ dl)cj){r) 

72(f) = d x d y 4>{r) (14) 


The convergence field takes the form 


«(f) = ^ V 2 </>(r) = i(<9 2 + dl)(f>{r) 


( 15 ) 


where V 2 = did 1 = didi as the metric is given by g = diag( 1,1). 
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C. Weak lensing shear as a spin weight 2 object 


Weak lensing induces a complex shear field 7 (r) which transforms under a rotation (by an angle t/j in the anti¬ 
clockwise direction in our conventions) of the fixed coordinate system according to 7 —> 7 e _ls ^, where s = 2 is its 
spin weight. This phase dependence expresses the fact that the complex shear field 7 ( 7 ") is invariant under a rotation 
over 7 r radians. Also, the two components of the shear field, 71 and 72 , are related by a 7 t /4 radians rotation which 
transforms one field into the other: = —72 and 72 = 71 where the prime denotes the transformed fields. In this 

subsection we use a description of the field based on a geometrical differential operator called edth and symbolised 
by 9. This operator was firstly introduced by Newman and Penrose in 1966 [1(1 (see also |56|) on the surface of the 
sphere in order to define the now widely used ‘spin-weight s spherical harmonics’ which live on the 2D spherical space. 
A spin-weight s spherical harmonic, symbolised as s hra, can be seen as a generalisation of the scalar, vector and 
tensor spherical harmonics (see [Ell 03)- The s hm (not defined for |s| > l) form a complete orthonormal basis for 
each s, a result that we will use later in order to describe the spin-weight s functions on the 2D sphere. The spin s of a 
function is related to its transformation properties under a rotation of the frame field where it is defined (for example 
a scalar is a spin-0 function). The 9 (and its complex conjugate 9) act as raising (and lowering) operators on the 
’quantum number’ s (integral number), such that s Y) m are eigenfunctions of 9 9 and can be obtained by applying 
9 to the standard spin-weight 0 spherical harmonics Y) m . The operator 9 is effectively a covariant differentiation 
operator acting on the surface of the sphere which enables one to relate quantities of different spin. In particular, 
it allows one to conveniently relate spin-weight s objects, which are not invariant under rotations of the coordinate 
frame, to scalar quantities, which are invariant quantities under rotations. All invariant differential operators on the 
sphere may be expressed in terms of it. 

The operator 9 and the concept of spin-weight s functions have been generalised to any 2D Riemannian manifold 
(a 2D manifold with a metric) (e.g. |~fT3. For details and further references see the Appendix, in particular 

Eq. jAl for the expressions relating the 9 operator to covariant derivatives in a 2D Riemannian manifold. 

As pointed out originally by Newman & Penrose pf| . any spin-weight s function defined on a 2D Riemannian 
manifold can be uniquely decomposed into a scalar gradient (or electric/even) ‘E’-component and a scalar curl (or 
magnetic/odd) ‘B’-component. We say that the spin-weight s field 77 is even if 77 = 9 s / and odd if 77 = i 9 s / for 
some real-valued spin-weight 0 function / (if s < 0 then we interpret 9 s as (— 9)I s !). For s = 1 this decomposition 
corresponds exactly to the classical Hodge-Helmholtz decomposition of a vector field into the sum of a gradient and 
a curl component. 

As a spin-weight 2 object, the shear field can be written as the second edth-derivative of a complex potential 25, 

ssaiSssiii: 

7 ( r ) = 7 i(r) + * 72 (r) = ^9 d[c/> E (r) +i(j> B {r)\, (16) 

7*(r) = 7i( r ) - *72 (r) = ^9 9 [<Mr) - *<M*0] (17) 

where we introduce two scalar real functions 4> E ( r ) and 4>b (r) for the even and odd parts. The normalisation factor 
of 2 was chosen so that one can later identify immediately the lensing potential to the even field (f> E ( r )- As we will see, 
weak shear is derivable from a real (lensing) potential <f> E (r) = requiring $3(r) = 0. This definition differs from 
the convention used in CMB polarisation studies (see for instance j43|,[44|) and comes about due to the conventions 
adopted for the spin-weight formalism (see the Annendix IA II) . Note that the scalar functions (j) E and <f>B introduced 
are invariant under rotations of the reference frame. 

We emphasise that the spin-raising and lowering operators 9 and 9 act on the 2D manifold at a distance r. We are 
interested in their expressions both in the full-sky 2D sphere and in the flat-sky 2D Euclidean space. Their derivation 
in both geometries for any value of spin s is detailed in the Appendices IA 21 and I A 31 and we give the final expressions 
for 9 and 9 here. In the 2D spherical full-sky we have 


S sf(9,<p) 


and 

9 sf(0,<p) 

In the 2D flat-sky the latter reduce to 

9 

9 


= — (dg + icscOdcp — scot 6) s f(9, <p) 

= — sin 5 9 (dg + icsc 9 d v ) sin -s 9 s f(9, y>), 

(18) 

= — (dg — icsc9d v + scot 9) s /(0, ip) 

= —sin _s 9 (dg — icsc 9 d v ) sin s 9 s f(9 , ip). 

(19) 

sf(x,y) = ~(d x +id y ) s f(x,y), 

a f(x,y) = -(d x -id y ) s f(x,y). 

(20) 
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If we use these results, or rather these results expressed in terms of covariant derivatives, Eq. m can be re-expressed 
explicitly in terms of covariant derivatives on the sphere and in the tensorial formalism used in section htbi as 


[7 ( r ) ]ij 


7i (r) 
sin 072(7") 


sin 072(7") 
—sin 2 071 (r) 


\ [V e V 0 - csc 2 0V v V„] 

—esc 0 V ^ V e 
i[sin0 VeVg — csc0 V<+ 


v„v e 

[V v V v - sin 2 0 V e Ve] 

i[sin0 VgVe — csc0 V 
1 sin 0 XI ..XI a 


<M r 




( 21 ) 

<t>B{r) 


where we obeyed to the conventions defined in the Appendices EH and El In particular, we projected the complex 
cosmic shear 7 into a tensor [7 ] ^ by means of relation Eq. El and again used the general result VgV^ (j> = V^Vg </>. 
We can now compare directly the general E-B decomposition of the shear Eq. (ZD with the theoretical predictions 
for gravitational weak lensing shear on the 3D full-sky Eq ©■ We can do this because the choice of the 2D polar 
coordinate system in both cases is the same. It is then mathematically straightforward to identify 4>e with the lensing 
potential <f> such that 4>e( 7 ") = (j>(r) and </>_ b ( t ") = 0 . Similarly one can make such an identification in the flat-sky 
approximation and re-express explicitly Eq. m in terms of partial derivatives on a 2D Euclidean space in the same 
tensorial formalism as above 


[7( r )k = (didj ~^6ij\7 2 )4>E(r)+ ^{£kjdid k +£kidkdj)(l)B(r) (22) 

where £7 is the anti-symmetric 2-dimensional Levi-Civita tensor for which £00 = £11 = 0 and £01 = — £10 = 1 |26i |. If 
we compare this equation with the weak lensing theoretical expression in the flat-sky limit Eq. d, the identification 
between <f> and <f>E is, as expected, immediate. Therefore the shear field induced by grav itational tidal fields only 
produces an E-pattern in cosmic shear maps. Such a result has been shown before (H, HI, H3, IM, |65| and has been 
applied to weak lensing data in the flat-sky limit 0, A, (| 0, @ ■ It was to be expected as density (scalar) perturbations 
only produce E-type effects. Hence the decomposition of the weak lensing data into E and B components presents a 
clear advantage over the 71 and 72 decomposition. On the full-sky it allows one to isolate the effect caused by weak 
lensing and disentangle it from any non-lensing contributions like noise, foregrounds or systematics which should 
contribute similarly to both the E and the B modes. It is worth noting that both gravitational waves [27), multiple 
light lensing scattering effects and source clustering 0] may also induce B modes, but their levels are expected to 
be small. As we will mention shortly, one can have additional spurious curl-type effects caused by the observational 
pixelisation or by finite fields, due to a leakage between the E and the B modes [43j, j+j]. In such a situation the 
B-rnode acts as a vital test for any non-gravitational signal in the data. Aside from such complications arising by 
finite fields and pixelisation, the E and B potentials have also the advantage of being rotationally invariant and so 
there are no ambiguities in their definition related to the rotation of the coordinate system on the sphere. Their 
interpretation is thus clearer than the one of the 71 and 72 decomposition over the whole sky. 

Having established this correspondence in 3D, one can recast both Eq. 0 and Eq. d. as the following valuable 
relation between the weak lensing shear 7 and the lensing potential in terms of the 8 and 3 operators 


7 (r) = ^3 3^(r), (23) 

7 *(r) = ^dd<P(r) (24) 


or, alternatively, the equivalent relation for the two orthogonal components of 7 

7i( r ) = ^(3 3 + 3 3 )</>(t"), 

72 (r) = -^(3 3-33 )<j>{r). (25) 

Similar expressions can be obtained for the scalar convergence field which is a spin-weight 0 object. As shown 
in the Eq. (IA8II of the Appendix, the Laplacian acting on a rank-0 tensor is equivalent to applying the combination 
[3 3 + 3 8]/2 to the corresponding spin-weight 0 quantity. In particular, if one compares directly Eq. IA1 41 expressed 
in the full-sky to Eq. o one sees that Eq. m can equivalently be written as 


n(r) 


-[8 8 + 3 3 )4>{r). 


(26) 


Such an identification can also be performed in the flat-sky case. Note that there are natural extensions to higher 
spin-weight objects: E(r) = —^(338 + 333 + 833 )(j>{r) and Q(r) = —58 8 are the analogues of the first 

and second flexion respectively (see [gTL f6^| j. 


We now have both the weak lensing shear and convergence fields expressed in the edth formalism which encompasses 
both the full-sky and the flat-sky limits. Due to the role of the 3 (3 ) differential operators as spin-raising (lowering) 
the spin of spin-weight s spherical harmonics, this alternative compact mathematical tool will be of particular beneficial 
use when studying weak lensing in the 3D full-sky. 


III. DECOMPOSITION OF 3D FULL-SKY WEAK LENSING INTO SPIN-WEIGHT SPHERICAL 

HARMONICS 

In the previous sections we have expressed the weak lensing shear and convergence fields as spin-weight functions. 
Aided by the edth formalism we have established their relation with the weak lensing potential by means of a simple 
and unique E-B decomposition scheme. In the present section, we concentrate on the 3D (geometrically flat) spherical 
sky where the previous fields can be decomposed into a combination of spin-weight spherical harmonics and Bessel 
functions. We thereby derive the relation between the expansion coefficients of the weak lensing quantities and the 
expansion coefficients of the weak lensing potential. The latter is readily related to the expansion coefficients of the 
gravitational potential and hence to those of the matter density field. 


A. Representation of a general distortion field in spin-weight spherical harmonics 

For spin-weight 2 fields such as the weak shear field, the natural 3D basis functions are products of radial functions 
and spin-weight 2 spherical harmonics. With this choice, the 3D expansion coefficients of 7 are related very simply to 
the expansion coefficients of the lensing potential, where the latter is expanded in products of ordinary (spin weight 
0) spherical harmonics and the same radial functions. 

Let us first assume a general distortion shear field rj(r). This field can be decomposed in two orthogonal components: 
r)(r) = rj\ (r) +irj 2 (r). Because of its transformation properties, such a field is a spin-weight 2 object. By consequent 
it can be decomposed into an even/‘E’ and an odd/‘B’ part by means of two scalar real functions defined on the 2D 
sphere of radius r, </># (for the even part) and </>b (for the odd part) as in Eqs. (FTTH) and 

V{r) = ^3 3[<Mr) +i<t>B{r)] , rf{r) = ^ 3 3 [0 E (r) - i<t>s(r)}. (27) 

As discussed previously, the two real scalar functions 4>e and (j>B introduced completely characterise the distortion 
field. They have the advantage of being scalars which are invariant under any rotation of the coordinate system. We 
can then expand them in 3D on the full sky as we expanded the lensing potential field in Eq. in terms of the 
standard (spin-0) spherical harmonics Y( m and of Bessel functions 


r°o 00 ' 1 ' t£ _ 2 )! 

<t>E{r) = -2 / dk^2 \ (o 1 9(1 a E,em{k) Z ktm (r, 9, ip), 

<=0m=-<V ^ + 

<t>B{r) = -2 [ dk^2 ^2 \ a B,im{k ) Z ktm {r, 0, y) (28) 

<=0m=-iV ^ + 

where we introduced the orthonormal basis functions Z k ( m for simplicity and which, in a spatially flat background 
geometry, are given by 


Zum(r,6,ip) 



kj e (kr)Y em (e,ip). 


(29) 


Apart from the normalisation factor of 2, the normalisation terms and the sign in Eqs. (EH were introduced to maintain 
consistency with previous works on CMB polarisation |44|. The factor 2 is linked to the choice of the spin-weight 
conventions as established in the Appendix^and is necessary if one wants to recover relations between the expansion 
coefficients of 7 and of the fields <pE and (j>B which are analogous to the relations between the CMB polarisation field 
Vij and the electric and magnetic scalar fields |44|. We point out here that we do not assume 3D homogeneity and 
isotropy (at a fixed instant in time) for the scalar fields 4>e and cj>B and so do not introduce the time r dependence in 
the coefficients of their expansion as discussed in section EU If we now substitute Eqs. E3 into Eqs. am and use 
the properties of the operators 3 and 3 acting on the standard (spin-weight 0) spherical harmonics (see Appendix 
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Eqs. 11661) 1 one obtains naturally the following relations for 77 

r°° 00 1 

v( r ) = dk EE 2 rUm{k) 2 Zum{r, 9, ip), 

"'° e=0 m=-l 

pOO 00 ^ 

7 j*(r) = dk EE -JHm{k)- 2 Z k tm(r, 9, ip) (30) 

, ' 0 t= 0 m=-e 

where the expansion coefficients are related to the <j>E and <ps expansion coefficients by ITU I3| 

2 V^m{k) = [&E,£m T ^B,£m\ (k) , _‘/]£m(^k) = [0'E 1 £m ^B,£m\{k) j (31) 

1 X 

— 2 [ 2 V£m H - —•^i£m\{}^) 1 — 7) [ ^f]£m (&)• (32) 

and where now the orthonormal basis functions ± 2 Zk£m are expressed in terms of a set of functions, the spin-weight 
±2 spherical harmonics ±2 Yim, the ideal basis to express any spin-weight ±2 function 

kje(kr) s Ye m (0,ip) (33) 

where the spin is given by s = ±2. As we can see, any distortion field 77 defined on the 3D full-sky is most naturally 
expressed in terms of spin-weight 2 spherical harmonics, which define a set of orthonormal basis on the surface of the 
sky. The edth formalism introduced in the previous section has just shown its mathematical advantages when working 
in spherical space with spherical harmonics. 

Inserting Eqs. ED into Eqs. CT we can also relate the 771 and 772 orthogonal components to the E-B decomposition 
expansion coefficients of the distortion field as 


,Zu m {r,6,ip) = \ — 
7 r 


Vi(r) 


m{r) 



00 £ 

dk EE [aE,tm(k) Xi ,k£m (r,9,ip) + ia B ,em(k)X 2 ,k£m (r, 9, ¥>)], 

£=0 m=—£ 



00 £ 

dk EE fa B ,£m(fo) X 1 

,k£m (r, 9, ip) 

£=0 m=—£ 


iaE,£m(k)X 2 

,k£m (r, 9, <,?)] 


(34) 


where we defined two new basis sets: X ltkim = ( 2 Zum + - iZum)/ 2 and X 2 ,k£m = (2 Z k £ m - _ 2 2 Mm )/2. As we have 
the relations X^ Mm = -Xi^-im, x 2Mm = -A \k-im, a* E tm = a E -im and a* B tm = a B -im, the coefficients 71 and 
72 are real quantities. 

Both representations can be chosen to study a general observed distortion field 77, i.e. either use the 771 and 772 
components or transform them into E and B type quantities. They contain the same information but, as explained 
before, the E and B are rotationally invariant on the full-sky and behave differently under parity transformation 
ITa lia . In terms of the expansion coefficients, a,E/m(k) has parity (—1)^ while a B ,i m {k) has parity (—l) f+1 , a 
property which may certainly be useful when trying to characterise the origin of the distortion field observed on the 
sky. In particular, the E-B decomposition components have a special meaning in the context of weak lensing studies as 
weak lensing is in nature an E-type field and does not produce any B-type component, as we saw before in section If T Cl 
In practice, if one wants to determine the real scalar functions (f>E and cf>B describing a given distortion field on 
the sky one should use Eqs. ED- This is a trivial task if our observations are performed over the full-sky because 
one can use the orthogonality of the spin-weight spherical harmonics to obtain the E and B expansion coefficients 
from the measured field 7 . But full-sky observations are unrealistic. Even if an experiment does cover the whole 
sky, certain regions always need to be removed in order to minimise the foreground contribution (such as bright 
stars). The perverse side effect of partial sky-coverage is to induce an E-B mode mixing. In this case the non-local 
decomposition <ED is not unique. If one tries to solve for 4>e or <j>B by taking linear combinations of second derivatives 
of the distortion field one gets 

V 2 (V 2 + 2)(j)E = ^[9 977+ 3 3’?*], 

V 2 (V 2 + 2)0 b = ^[33,7- SSt?*]. (35) 

This result enables one to use distortion observations to construct the E and B distortion shear modes on the full-sky. 
But over a cut-sky, solving this system implies the specification of unknown boundary conditions (i.e. value of the 
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potentials and their derivatives at the boundaries). This problem has been extensively studied by many in the context 
of CMB polarisation measurements UJ] but to date separation of E-B modes has been performed statistically 

during the estimation of the power spectra using the so-called pseudo-Q method (see @). On the flat-sky and 
in the context of weak lensing, the locally defined aperture mass M ap may similarly provide unambiguous mode 
separation [67ll6^ |. A similar statistic can be built for the weak lensing B mode. It is beyond the scope of this work 
to enter into technical details related to mode separation,so we refer the reader to Challinor j^| for a useful recent 
CMB polarisation review and van Waerbeke and Mellier |2| for an updated weak lensing review. 


B. Representation of the weak lensing shear and convergence fields in spin-weight spherical harmonics 

The previous results can easily be applied to the case of the weak lensing shear field 7 for which (see Eqs. lt23l) 
and J0J) 

7 (r) = igg^(r) , 7* (r) = ^9 d<t>(r) 

where <f) is the lensing potential defined in Eq. 0 and which we can expand in 3D as in Eq. 0 where the expansion 
coefficients are <fie m (k ). Comparison of the 3D expansion of the lensing potential (f> with Eq. (1281) shows that the a E /m 
and the <f>e, m are related by 


a E ,tm{k) = (g _ 2)i ^ m (^)~ ( 36 ) 

As demonstrated before, for gravitational weak lensing, a,B,em{k) = 0. The 3D full-sky coefficients of the spin-weight 
2 shear field 7 , 2 t e m (k) and _ 2 'ye m (k), defined similarly to Eq. (13(71) are then related to <\>t m by (see |24. I3& |4^| ) 


2 7 em(k) 


. 2 lim{k) 


1 l(e + 2 )! 

2 y (£ — 2)! 


fyzm (k). 


(37) 


The standard components of the weak lensing shear 71 and 72 can be directly obtained from this last relation. If one 
uses the Eq. © with aB 4 m{k) = 0 and the following expansion in 3D for the shear components 71 and 72 


r°° 00 *■ 

7i(r) = dk EE 

1=0 m=—i 
/»oo 00 ^ 

72 (r) = dk EE 

-'° t=0 m=-e 

where Xi^ m and X 2 ^kim were defined as in Eq. 13411 . one immediately has 

7i ,im{k) = ij 2 /m(k) = -a E ^ m (k), (39) 

where the expression for a E ,tm in terms of the lensing potential was defined in Eq. <29- As we see, aside from the 
physical and experimental advantages of the E — B decomposition scheme, using either of the expansion coefficients 
7i/2,tm or flBtm is equivalent as they are very simply related. 

Likewise, we can do the same exercise for the convergence field n. From Eq. 03) and using Eqs. IB(il) . the (scalar) 
expansion coefficients of k are given by |t^| 

«/m(*) = -E±ll^ m (fc). ( 40 ) 

To summarise, the coefficients of the expansion of the shear field in spin-weight 2 spherical harmonics and spherical 
Bessel functions, ±2 7^ m (fc) are related to those of the lensing potential 4>i m (k) by Eq. ifTTI) . The (j)e m {k ) are related 
to the gravitational potential <E , f m (fc) by Eq. 0, and these are in turn related to the overdensity field by Poisson’s 
equation ©. In this way we establish a connection between the observable shear quantities ±2 7 e m (k) (and associated 
quantities) and the overdensity field, whose statistical quantities have known dependence on cosmological parameters. 
In the knowledge of such results one may proceed to the derivation the 3D full-sky power spectra (or any any relevant 
statistics) of the weak lensing shear and convergence fields. 
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IV. WEAK LENSING 3D POWER SPECTRA 

We shall now be interested in calculating the weak lensing shear 7 and convergence k 3D power spectra and in 
relating them to the 3D gravitational potential power spectra both in the full-sky and flat-sky approximation. 


A. Weak lensing 3D power spectra on the full-sky 

If one performs a 3D spectral decomposition in the full-sky of a statistically homogeneous and isotropic field /(r; r) 
at time instant defined by r (one can use the comoving distance r as a measure of real time t as r and t are physically 
equivalent being related via the scale factor a) with coefficients fe m {k',r ) as in Eq. Q, then the 3D power spectrum 
of the field at time r, Q(fc;r), is defined by 

(, fem(k;r)f£, m ,(k';r )) = C e {k;r)5 D (k - k')6 U 'S m m' (41) 

where Sij is a Kronecker delta function and 5d(x) a ID Dirac delta function. The time dependence (introduced 
through r) of the coefficients, at which the spectral expansion is performed, was explained in section and is 
necessary to ensure the 3D spatial homogeneity (and isotropy) of the field we are expanding at a certain time instant. 
Most importantly, Ce(k; r ) is independent of £, and is simply the 3D power spectrum P(k ; r): 

C e (k-r) = P(h, r) (42) 

where P(k;r) is obtained from performing a 3D Fourier expansion of the held /(r;r) so that (/(fc; r)/*(fc , ; r)) = 
(2ir) 3 P(k]r)5 3 D (k — k'). Indeed the spherical expansion in 3D of Eq. © for f(r;r) is equivalent to the 3D Fourier 
expansion 

f(r;r) = ^3 J d 3 kf(k; r) e lk r (43) 

because of the well-known identity j4fj 


m=£ 

e ik.r = 47 t^ ^2 i e jt(kr)Ye m (0 k ,ip k )Ye m (9 r ,ip r ). 

£ m=—£ 


(44) 


One can see this by substituting Eq. (EU into Eq. m and using the identity d 3 k = k 2 d£l k dk where f 1 k = (9 k .tp k ) so 
that we obtain the following relation between the expansion coefficients fim{k\r) (see Eq. @) and /(fc;r) 


hm{k\r) 



dQ k f(k;r)Y em (9 k ,ip k ). 


One can then substitute this last expression in Eq. 1411) and use (/(fc; r)f*(k'\r)) = (2n) 3 P(k\ r)8 3 D {k — k'). The trick 

then relies in replacing 5 3 D (k — k') = J d 3 ^' x /( 2tt ) 3 and using the identity 1441) . The orthogonality relations 
of the spherical harmonics 1B2I) and of the Bessel functions cost as well as Y^ m = (—1 ) m W-m will then be useful to 
reach the final result, i.e. Ct(k\r) = P{k\r). 

Now if we use the expansion coefficients of the field located at different distances (or equivalently at different times 
from us) say at r and r', to effectively calculate the cross-power spectrum of two different homogeneous and isotropic 
fields, /(r;r) and /(r';r'), then one can still argue that the 3D homogeneity and isotropy argument holds: 


(a^fcjr)^,™,^';/)) = C e (k;r,r')5 D (k - k')5u'5 m m> 


(45) 


where now we have the identity 


Ct{k\ r, r') = P(k) r, r') (46) 

where again P(k;r,r') is defined by (f(k;r)f*(k';r')) = (2ir) 3 P(k;r,r')S^,(k — k'). A similar result is frequently 
utilised to determine the cross-power spectra of different homogeneous and isotropic fields, like for instance the weak 
gravitational lensing of the CMB and the cosmic shear Helds |53 • It will be necessary to us when calculating the power 
spectra of the gravitational potential field $ which comes about when calculating the spectra of the lensing potential. 
Again we stress that, contrary to the gravitational potential, the 3D lensing potential <f> is not homogeneous and 
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isotropic in 3D space. It is given by a 2D projection at each source distance r of the gravitational potential existing 
between us and the source and so it maintains the homogeneity and isotropy characteristics of the gravitational 
potential field on the 2D sky, but not in the radial direction. We will then use for the field <fi the relation 

= Cf(k, k')dwS mrn ', (47) 

where Cf^ is the 3D lensing potential power spectrum. Contrary to the case of a 3D homogeneous and isotropic 
field, an equivalent relation to Eq. 02 does not hold any longer. Naturally, the lensing potential, the shear and the 
convergence all share the same statistical properties. 

For completeness, although redundant, we present the various possible shear power spectra depending on the type of 
shear decomposition chosen. If one uses the expansion coefficients of the shear j(r) itself, i.e. ±2 'ye m (k;r) of Eq. (liTTl) . 
then one gets 

CT(kiM) = \ c t\k u k 2 ) (48) 

where Cf ^ is the lensing potential power spectra. All the remaining power spectra can then be conveniently expressed 
in terms of Cj 1 (k\, k 2 ). For the E-B decomposition components, i.e. a_B^ m (fc) of Eq. ll36l) and aB,tm(k), one has ill] 

Cf B (fci,* 2 ) = Cj\k u k 2 ), 

Cf B {kiM) = 0. (49) 

And finally for the 71 and 72 shear decomposition Eqs. 02 and 02 one obtains 

c7 i7i (fc 1 , k 2 ) = c7 272 (/c!, k 2 ) = ap{kx, k 2 ), ( 50 ) 


where for symmetry reasons C 7lT2 = 0. Note that the expansions of 71 and 72 are in terms of the orthogonal but 
not orthonormal functions Xi^ m and X-> fe / m , so 71 and 72 each contribute half to the power of 7 , as expected. 
Concerning the convergence field we have 14211 


cr(h,k 2 ) 


£ 2 {£ + l) 2 
4 


C^(Aq,fc 2 ). 


(51) 


The next and final step is to relate the lensing and the gravitational potential power spectra, Cf^ and 
respectively, by means of Eq. ©■ After some straightforward algebra we obtain 


cr&iM) = 


16 


k 2 dk h(ki,k)Ie(k 2 , k) 


(52) 


with 


I e (ki,k)=ki J drrje(kir) J dr' J j t {kr')yj P^{k\r') 


(53) 


where we have used Eq. 02 applied to the gravitational potential field and introduced the familiar 3D power spectrum 
of the gravitational potential via E( i. 02 - For a non-uniform clustering of sources, the number density needs to 


be introduced into the outer integral The correlations in the potential field are significantly non-zero for small 
separations (much smaller than the speed of light times the timescale over which it evolves), so we have assumed for 
convenience that the gravitational potential power spectrum can be accurately approximated by 


}<[><£ 


(fc; r, r') ~ y]p**{k\r) P**{k;r'). 


This can be seen by noting that Poisson’s equation implies that 

, 3 n m H 2 

$(fc; r) = - . “ 1 5(fe; r ) 


and so 




3n m H 2 


2 a(t)k 2 


d 3 rd 3 r' (5{k;r)6*(k';r')) _jfc. r+i fc' r ' 
a(t)a(t') k 2 k' 2 


(54) 


(55) 


(56) 
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and the correlation of S is restricted to small scales | r — r' |< lOOMpc. The lookback time over such a distance is 
small, so we can approximate r ~ r' (or t ~ t'). Therefore we can replace the power spectra P $$ (fc;r, r') by either 
r) or r'). For algebraic convenience, we choose the geometric mean of the power spectra, which allows 

us to separate two internal integrals and reduce computation time significantly. A further justification is apparent in 
Fig. ©: the Bessel functions cut off long-wavelength contributions with k < llr max where r ma x is the extent of the 
survey. Practical issues of implementation are briefly discussed at the end of the next section. Note that the two-point 
statistics of the shear field are dependent on the nonlinear potential power spectrum, which is related to the nonlinear 
matter power spectrum via Poisson’s equation. Accurate fits for the nonlinear matter power spectrum are available 
to enough wavenumbers |7fl , but the accuracy for high-precision determination of dark energy properties needs to be 
tested with simulations. 


B. Weak lensing 3D power spectra in the flat-sky limit 


In the flat-sky, we can expand a 3D field / at 3D position r = (r, 9) on the sky into a combination of 2D Fourier 
modes and Bessel functions in the radial direction r 

_ P? r™ r°° rP-f - 

/(A 0 ) = \ ~ / kdk / f(k4)ji{kr)e lL0 , 

V 7T 7o J o (27 T) 2 

/ i) /• oo r oo 

- / r 2 dr / d 2 9f[r,9)kjt{kr)e~ lLe 

n Jo Jo 

where the notation x describes 2D vectors. Such an expansion was chosen to maintain a direct relation with the 3D 
full-sky expansion introduced previously. It thus presents the same attractive advantages as its full-sky counterpart 
(see discussion in section nm . One can now easily obtain the full-sky to flat-sky correspondence, i.e. establish a 
relation between the expansion coefficients f(k,£) in the flat-sky as defined above and ft m (k) in the full-sky defined 
in Eq. iflf. 

For small angles around the pole of spherical coordinates, defined by the an gles (6, p) ignoring the radial dependence, 
for which 6 —> 0, one can use the following expansion of the plane 2D wave [42. l7ll| 


(57) 

(58) 


(59) 


where t = (£ cos pt, ism pf ) and 9 = (9 cos<p,9 svo.ip). The vector £ can be interpreted as the continuous limit of 
the integer that labels the spherical harmonics Yi m . Replacing Eq. m into Eq. m and making use of / 0 °° d?i = 

idl j‘j~ dpt. P2 1 £ / 0 27r dpt for high l, one obtains the following relation between the 3D flat-sky and the 3D 

full-sky coefficients 

hm{k) = ^i m j 2 J ^e- im *'f(k, £). (60) 

The inverse relation is 

/M = (si) 

m 

We point out that, as expected, we recover the same relations that were obtained by others for the 2D case EHInl , as 
the small-angle limit only affects the 2D angular expansion, not the radial direction. As the reader may have guessed, 
we also recover that the flat-sky and the full-sky power spectra are equivalent in the high~£ limit. For a 
fully 3D homogeneous and isotropic field, like the gravitational potential we have 

C^ikJ^r) ^=°° Cf®(k;r) =P $4> (fc;r) (62) 

where we have used Eq. 1411) . the definition ($(fc, /)$*(fc', £')) = (27r) 2 C $<I> (fc, £)So{k — k')5p(i£ — l') and Eq. 11421) . 
Again p $4> is the 3D gravitational potential power spectrum. For a 3D field which is homogeneous and isotropic only 
on the 2D sky (e.g. the lensing potential) 


C+*(ki,k 2 ,l) £ =°° Cf*(fci,fc 2 ) 


(63) 
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where in this case we have defined the power spectra via £)<f>* (k', £')) = (27r) 2 C^(fc, k', £)5 2 D (£ — £') and 

= Cf+ik, k')5 U '5mm'- 

One then has the following identical expressions for the flat-sky power spectra of the shear and the convergence 
fields ;42] 


C^(h,k 2 ,£) = (64) 

C™(h,k 2 ,£) = ^C^(k u k 2 ,£). (65) 

where now 

16 r°° 

C^(fci, ki,i) = - 5 -j / k 2 dkl(k i, k)I(k 2 , k) ( 66 ) 

tt 2 c 4 Jo 

with I as before (see Eq. 15111 1 1. The equations 1611) . 1^7)1) and open up the possibility of cosmological parameter 
estimation with the shear field, whose 2 -point statistics are here related to those of the underlying gravitational 
potential field (and hence to the matter density via Poisson’s equation). For clarity, we have not considered here the 
practical issues of noise, source selection function, weighting schemes, or errors in the photometric distance indicators. 
These can be dealt with, and are studied in a companion paper j&jj ]. which addresses the accuracy with which 
cosmological parameters can be estimated with current and future 3D weak lensing surveys. Also considered in that 
paper is that one has to assume a fiducial cosmological model in order to translate source redshifts into comoving 
distance coordinates r. 

In Figures 0 to 0 we illustrate the amplitude and shape of the signal produced by the 3D shear power spectrum 
C 77 (fci, k 2 , £) of Eq. fim one may ideally expect to measure. We also show how the shear power spectrum varies 
with the value of the equation of state parameter w. The plots are for a fiducial ACDM cosmological model with 
parameters Da = 0.73, D m = 0.27, fit, = 0.27, H 0 = 71 km/s/Mpc and equation of state of dark energy parameter 
w = —1. Figs. 0 , 0 and 0 show the features of the 3D shear power spectrum by taking various cuts through the 
3D (fci, k 2 , £) space. We consider a survey to r max = 5000 Mpc, which replaces infinity as the upper limit of the radial 
integral in Eq. lt53l) . The sampling interval in wavenumber k is 2.5 x 10 4 Mpc 1 and in multipole log w (£) is 0.025. 
We also illustrate in Fig. 0 how the shear power spectrum varies with w. The differences are slight in an individual 
power spectrum, but as this is a 3D study, there are many useful £ modes, which increase the sensitivity of the 3D 
shear as a cosmological parameter diagnostic. 


V. 3D CORRELATION FUNCTION 

For surveys with complicated geometry, the 3D shear correlation function may be a more appropriate tool to use, 
in place of the 3D shear power spectrum. It is straightforward to relate the 3D shear correlation function to the 
lensing potential power spectrum (and hence the matter power spectrum), via (15211 . as follows. The shear expansion 
as defined in Section fill FI implies that 

( 7 (r) 7 *(r'))= ^ f dk dk'( 2 j Mm 2 ^m(r, 9, tp) 2 Z* k , t , m , (r, 9', tp'). (67) 

Imi'm' 

Using 0U, this reduces to 

( 7 (r) 7 *(r')) = ^ Z (d — ^)! / dk' 0 ^{k, k')kk'j e (kr)j e (k'r') 2 Y em (d, ip) 2 Y e * m (6',ip'). ( 68 ) 

Using the generalised addition theorem for spin-s spherical harmonics l7Gl |, we find 

( 7 (r) 7 *(r')) = ^ 4 ^ 1 (^_2)! / dkdk ' C t <>> ( fc ’ k')kk'j e (kr)j i {k'r') 2 Y e _ 2 (J3, a)e~ 2ie , (69) 

where (3 is the angle between the directions (0, ip) and (O' , ip'), and a and e are the angles between the line joining the 
two points on the celestial sphere and the lines of constant ip. An explanatory diagram can be found in 0 . Note 
that one might expect that the correlation function would depend only on the angular separation through /?, but as 7 
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FIG. 1: The left-hand panel shows a grayscale cross-section through the C 77 (fci,fc 2 ,£) 3D matrix where each element has been 
multiplied by k\k 2 . The black/white colour represents positive/negative values. The right-hand panel represents a cut in the 
(fci,fc 2 ) plane for a given value of l and k\. The fiducial cosmological model parameters chosen are Da = 0.73, = 0.27, 

f lb = 0.27, Ho = 71 km/s/Mpc and w = —1. 



k/Mpc 


k/Mpc 


FIG. 2: These plots show diagonal cuts in the (fci,fc 2 ) plane at different £ values through the C 77 (fci,fc 2 ,^) 3D matrix. The 
effect of the approximate Bessel function inequality, kr > £, in Eq. 1531 can be seen. As the L value increases the diagonal 
terms of the covariance do not become significant until kr m ax ~ £, where r ma;E is the upper limit imposed on the r integral. 
The fiducial cosmological model chosen is the same as in Fig. (1). 
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t mode 


FIG. 3: Dependence of the diagonal part of the 3D shear power spectrum C 77 (fci,fc 2 ,^) on £, for various k values in Mpc 1 
The fiducial cosmological model chosen is the same as in Fig. (1). 


is defined with reference to the orthogonal coordinate system ( 6 , ip), there is an additional dependence on the angles 
a and e. Note also that the 3D correlation function is not homogeneous. If one chooses, one may simplify slightly 
by considering pair correlations of the shear defined with respect to axes parallel and perpendicular to a great circle 
joining the two galaxies on the sky. The rotation of the 7 coefficients (to 7 ) then gives 

{r')) = ^ 1 (g - 9 )i / dkdk ' c t^( k ^ k ') kk 'h(kr)ji(k'r')2Y e -2(P,0) (70) 

where 2^,-2 ($,0) may be simplified to a weighted sum of 3 associated Legendre functions. 

One can derive equivalent flat-sky relations by using the expressions m for the weak-lensing shear on the flat-sky 
and the expansion <E7|) for the lensing potential. Using the result 7 (k,t) = [(£ 2 — t 2 )/2 + i£ x £ y ] <f)(k,£), where £ and 
k are defined as in the previous section, we easily obtain 

(7M)7*(r',0')> = ^3 J d 2 ee 4 J dkdk'CM(k,tf,e)kVj i {kr)j t (k'r')e il -@- S> '> (71) 

where C^(k, k', £) is the flat-sky power spectra of the lensing potential. Alternatively, one can reach the same result by 
expanding Eq. around the pole of spherical coordinates 6 ~ 0 and taking the continuous limit of high multipoles 

e. 


VI. CONCLUSIONS 

In this paper we have developed in detail a complete formal study of 3D weak lensing, where one has estimates of 
the weak lensing shear field at known locations in 3D space. Most cosmological weak lensing surveys to date have 
been analysed in projection on the 2D sky, where the individual distances of the lensed galaxies are ignored. With 
distance information, a more sophisticated analysis is possible, which offers the prospect of greater statistical power. 
A partly 3D approach is to divide the sources into redshift slices, a process often referred to as tomography, but 
this crude division fails to explore the complete potential provided by distance information, and it makes sense to 
exploit as fully as possible the 3D information available: one essentially has a noisy estimate of the 3D shear field at 
the positions of all the source galaxies. For the estimation of some parameters, such as the amplitude of the power 
spectrum, 3D information adds relatively little, but for others, such as the equation of state of dark energy, 3D weak 
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k/Mpc k/Mpc 


FIG. 4: The derivative of the diagonal elements of the 3D covariance matrix C 77 (fci, for various £ values, with respect to 
the dark energy equation of state parameter w. The fiducial cosmological model chosen is the same as in Fig. (1) where w is 
allowed to vary. 


lensing analysis is very promising, and this facet of 3D lensing is explored in a companion paper j69| . Indeed, weak 
lensing on a cosmic scale may be the best cosmological method able to answer the important question of the nature 
of the dark energy. 

To effect a 3D analysis, we have used the theory of spin-weight functions successfully applied in the past to studies 
of the 2D CMB polarisation on the sky [S, Tijj. EH EH S3 > as well as theoretical 2D weak-lensing analysis [t| f26j, j27j. 
We have shown that the natural expansion basis to use for the 3D shear field is the product of spin-weight 2 spherical 
harmonics and spherical Bessel functions, which form a complete and orthonormal system of tensor functions on the 
three-dimensional space. Alternative expansions are clearly possible such as using tensor spherical harmonics j22| but 
these are more difficult to work with and have been in some way abandoned by the CMB polarisation community. 
With our basis choice, the two-point statistics (3D shear power spectrum, 3D correlation function) of the weak lensing 
3D shear field can be related in a straightforward manner to the 3D power spectrum of matter density fluctuations. 
The statistics of the 3D shear field are then connected to cosmological parameters via this relation, in addition to the 
redshift-distance dependence which enters the transform. This connection is easily shown using the edth differential 
operators. We are also able to relate the power spectrum of a full-sky 3D shear expansion to the power spectrum of a 
small-angle survey, where the use of spin-weight spherical harmonics is cumbersome, and a more familiar transverse 
Fourier expansion can be employed. This is the basis for cosmological parameter estimation using 3D weak lensing, 
and the formalism presented here provides a useful theoretical framework for 3D weak lensing studies. 
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APPENDIX A: SPIN-WEIGHT s FUNCTIONS AND SPIN RAISING AND LOWERING OPERATORS 

In this Appendix we review the notation and the mathematics of spin-weight s functions and of the geometrical 
spin raising and lowering (<3 and 9) operators defined over any two-dimensional Riemannian manifold (defined as a 
2D space with a metric). The relation between the spin-weight functions and tensor quantities is explained. We also 
clarify the relations between the 9 (9) operators and covariant derivatives on the 2-dimensional manifold. We then 
particularise to both the 2-dimensional unit sphere and unit Cartesian space, which will be of use to us in the main 
body of this work. We finally give an overview of the spherical spin-weigh t sp herical harmonics. Many articles and 
appendices in the literature have been dedicated to this topic (see |4j, [43, |45j, Ha, (5(J |23 to cite just a useful few). 

The spin-weight functions and the so-called edth (9 and 9) operators were first introduced in the 60s by Newman 
& Penrose j46j and further explored by Goldberg et al. |5(J. They mainly started to be used as a convenient tool in 
the context of the theory of gravitational wave radiation (see e.g. K.S. Thorne jjjjJ) an( i later they were introduced in 
the study of the Cosmic Microwave Background (CMB) all-sky polarisation as an alternative to the tensor formalism 
(see |4l|h The similarity between the CMB polarisation Stokes parameters Q and U and the weak lensing components 
7 i and 72 as well as the perspective of future large-field high resolution observations triggered recent full-sky weak 
lensing studies using the spin-weight formalism | 2 J, [42, . 


1. Representation of spin-weight functions in 2D Riemannian manifolds 

Let us define an orthonormal basis at any point of a two dimensional manifold with a metric gij {ei,e 2 }. The 
choice of the basis is not unique. Indeed, if we define the vectors m and m with respect to {ei,e 2 }, the choice of 
the basis is described up to a phase such that 

m = —=(e±+ie2)—^e~' l ^m, 
v 2 

m = —=(ei — ie 2 ) — * e l ^rn. (Al) 

y/2 

The vector m and its complex conjugate m obey the relations m J rrij = rn J rrij = 0, rrij = 1 and m 1 rrij + 
rn l rrij = <5*, where <5® is the Kronecker delta function. Here and throughout this manuscript, we use the Einstein 
summation convention. 

A complex function defined by s /(n) = fi(n) + i /2 (to), where f\ and /2 are real quantities defined on the manifold, 
is said to have spin-weight s if under the transformation Eq. (IA1I) it transforms as s /(to) —> e~ ls ^ s f(n). If we 
consider an arbitrary vector on the manifold v(fi) then, for instance, the quantities v.&i ± iv.e 2 transform as spin- 
weight ±1 quantities. Generalising to a rank-s tensor F il _j s (where i = 0,1), the object F ix ia m' 1 ... m‘” transforms 
as a spin-weight s object as each individual vector m ln contributes with a factor e _I ^ when it is transformed. The 
rank of the tensor is thus reflected in the transformation properties, and thus the spin, of the corresponding complex 
quantity constructed. 

One can therefore associate to every symmetric and trace-free component of a rank s tensor a spin-weight s quantity 
such that 


s /(n) = F il ... is m ll ...m ls , 

- s f(n) = F il „, is (A2) 

The trace-free condition refers to the vanishing under contraction of any two indices in the tensor. For example, for 
a rank-2 tensor, it refers to F\ = F it i g l 1 . As irreducible tensors of rank s in 2 dimensions have only two linearly 
independent components, say -F 00...0 an< ^ -P 10 ...O) we can combine them to create the associated spin-weight ±s complex 
quantities ± s f(n) = i 7 oo...o ± *J 7 io...o. It is a one-to-one mapping. Conversely, one can express any symmetric and 
trace-free rank s tensor in terms of spin-weight ±s objects 

FiL.A. = s f(n)rn il ...m lB + _ s /(n) m^... m ie . (A3) 

An important property of spin-weight objects is that there exists a geometrical spin-raising (lowering) operator 9 
(9) that have the ability of raising (lowering) the spin-weight s of an object such that under the transformation m 

9[./(n)] - e-^+^g [s /(n)], 

9[ s /(n)] - [ s f(n)}. 


(A4) 
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Most relevant is that these spin-raising and lowering operators are related to the covariant derivatives of the rank-s 
tensor associated to the spin-weight s object on the manifold 

3[ s /(n)] = -V2rn ll ...rn la rn k 'VkF il ...i s , 

3[ s /(n)] = -V2 m 11 ...rn l ‘rn k S7 k F h ,,, is (A5) 

for s > 0, and with m'“ replaced by m 1 " for s < 0. The choice of the sign and the normalisation are conventional but 
were chosen so that the form of the spin-raising and lowering operators on the 2D sphere match the original definition 
of Newman & Penrose [46} (see next section). The covariant derivative is defined as usual. For a rank-1 tensor Xj, 
for example, V jXi = djXi — T^Xk with the Christoffel symbol depending on the metric. If one replaces Eq. (IA3I) 
into Eq. ifATjli . one can rewrite Eq. (lA5l) for s > 0 as 

3[ s /(n)] = -V2 (X k [sf(n)]m k + st s f(n)) , 

3[ s /(n)] = -V2 (X k [sf(n)]rn k - st s f(n)) 

(A6) 

where 

r = V k m, m l m, k . (A7) 

The equivalent relations for s < 0 follow from these. One can express any differential operator in two dimensions in 
terms of 3 and 3 . For instance for a spin-weight s field s /, we have the following relation between the 3 operators 
and the Laplacian operator 


V 2 U/] = ^[3 3 + 3 3 } s f. 


(A8) 


2. Representation of spin-weight functions on the unit sphere 

We choose the orthonormal basis on the sphere to be aligned with the coordinate basis vectors 0 and (p of a spherical 
polar coordinate system {eg(n), e v (n)} where h is the radial direction vector normal to the surface of the sphere. 
The transformation defined in Eq. m corresponds to a right-handed rotation by an angle ^ of this basis around 
the vector fi. The convention chosen for the rotation is the same as in Zaldarriaga & Seljak ;41] and Okamoto & 
Hu [13] but differs from Newman & Penrose an d Goldberg et al. j5|| so that it conforms to the majority of CMB 
polarisation publications. 

The metric tensor is given by 


9ij ~ (0 sin 2 9 ) 

and the nonzero Christoffel symbols are 

r° n = — sinflcos# , = r} 0 = cot 6. 

The tensor components of the vector m (see Eq. ED) are chosen to be 


m 


= —— (1, icsc O') 

V2 


m k = — p(l, isin0). 
v 2 


(A9) 


(A10) 


(All) 


To obtain the exact relation for the spin-raising and lowering operators 3 and 3 on the sphere we use Eqs. (IA6I) 
and (1X71) . To calculate r we need the evaluate explicitly the covariant derivatives V,; rrij on the spherical basis with 
coordinates (#,<£)• We obtain 


V e m v 
X v mg 


0 , 

1 

—=cos 6, 

V2 

^^sin 6 cos 9 

V2 


(A12) 


such that 
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—cot# 

ET 


(A13) 


yielding Eqs. m and (frail . We note that the expressions obtained are identical to expressions found elsewhere HI 
|43, |45l |46}. This stems from the choice of a suitable formalism convention. Another useful result is the full expression 
of the Laplacian acting on a spin-0 quantity <f> 


l r - ^ 


= -[3 3 + 9 9]</> = [V 0 V 0 + esc 6* 


(A14) 


We remark that V 2 corresponds to — i(l + 1) in spherical harmonic space. 


3. Representation of spin-weight functions in Euclidean space 


We will also be concerned with fields defined over a 2-dimensional Euclidean space in the so-called flat-sky approx¬ 
imation for which the natural basis is the Cartesian coordinate system {e x (ft), e y (h)} where ft is the vector normal 
to the surface of the sky. In this case, the metric tensor is given by 


9ij — 



(A15) 


The tensor components of the vector m (see Eq. ED) are chosen to be 

= (A16) 

where m k = irtk- The Christoffel symbols are all zero and thus in Eqs. IA6|) and m we have that r = 0. In 
this case, the differential operators 9 and 9 naturally reduce to Eqs. m I 20 J. For completeness, the Laplacian in 
Cartesian coordinates reduces to 

v 2 0 = i[9 9 + 9 9]^, = [dl + d 2 y \t. (A17) 


APPENDIX B: SPIN WEIGHT SPHERICAL HARMONICS 


A scalar field defined on the sphere can be expanded in spherical harmonics Y( m (9, ip), which form a complete and 
orthonormal basis on the sphere. This basis of spherical harmonics is no longer appropriate to describe objects of 
spin-weight s. There exists similar sets of functions defined on the sphere that can be used instead. These are called 
spin-weight spherical harmonics s Ye m (9,p) and are defined by 0 


aYlm% if) = \j^^-D l _ s rn {9, ip, 0) (Bl) 

where D is the Wigner-D function (for technical details see Using the previous relation one can prove that 

the spin-weight s spherical harmonics satisfy the relation of orthogonality 




dcos 8 s 'Y ( *, m , (8, p) s Yim{8, ip) 




(B2) 


where 6ij is a Kronecker delta function. Also s Y ( * m = (—l) s _ s Y(_ m . A given spin-weight s function on the 3D sphere 
can thus be expanded as |49j 


sf(x) 


n s,£m(k) 



OO i 

dk EE sZktmipC') ®i 

£=0 m=—l 


/ 


d 3 x[ s f(x)} 


sZu m {x,8,ip). 


(B3) 
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where in spatially flat geometry the basis functions s Zktm are given by combinations of spherical Bessel functions 
and spin-weight s spherical harmonics as in Eq. (1351) . They are orthonormal 


J d 3 x s Z Um (x,6,y) s/ Z^ e , m ,(x,e,(p) = S D (k - k')S U 'Smm'S s 
where we have used Eq. G23 and the relation 


x 2 dx 


- kje(kx) 


- k' je(k'x) 


= 6 D {k~k'). 


(B4) 


(B5) 


If we are not assuming a spatially flat space then the radial functions need to be changed. For an open Universe for 
example je(x) —> X((Qk,x) where Xf is the hyper-spherical Bessel function (see Liddle & Lyth ( 4 ^ for details). 

A few properties of spin-weighted spherical harmonics when acted by spin-lowering and rising operators (as defined 
in Eqs. m and m) can be very useful 

3 s4{m = [{(■ ~ s)(£ + S + 1)] 2 s+ll^rra, 

3 s^lra = ~ [(^ + s)(i — s + 1)] 2 s-lYtmi 
3 3 s^tm = ~ ~~ s ){(- + s + 1 ) sYlmi 

3 3 sY tm = -(e+s)(e-s + l) s Yt m . (B6) 


In particular, the spin-0 and spin-2 spherical harmonics are related through 

a aw _ [Wji v 
3 3 Y(.m — W 2 )! 2 ^ m ’ 

3 3V _ / (^ + 2 )' w 

9 3 ^tm — W ^ _ 2^1 -2 Yim- 

A useful consequence of these properties is 


3 = ggggy, m = |i±|l|y, m = (£ + 2)(f + i)£(f-i)y fa 

which reduces to (see e.g. [44) j 

gggg = gggg =v 2 (v 2 + 2) 

when acting on spin-zero variables and where V is a covariant derivative on the sphere. 


(B7) 


(B8) 
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